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We investigate the low-energy quantum transport in ballistic normal metal-insulator - 
superconductor junction exposed to a magnetic field creating Abrikosov vortices in the supercon- 
ducting region. Within the Bogolubov- de Gennes theory we show that the presence of the subgap 
quasiparticle states localized within the vortex cores near the junction interface leads to the strong 
resonant enhancement of Andreev reflection probability, and the normal-to supercurrent conversion. 
The corresponding increase of charge conductance is determined by the distance from the vortex 
chain to the junction interface, which can be controlled by the applied magnetic field. 

PACS numbers: 74.25.QP, 74.25. Fy, 73.40.Gk 



I. INTRODUCTION 

Recently, the investigation of transport properties of 
normal metal- superconducting (N/S) hybrid structures 
has attracted a considerable interest. In the classical 
work by Blonder, Tinkham and Klapwijk-i it was shown, 
that at the energies below the superconducting gap Ao 
the charge transport can only be realized via the Andreev 
reflection at the N/S interface. This is a two-particle 
process in which the electrons with low energies e < A 
incident from the normal metal (N) are reflected at the 
N/S interface as holes traversing the backward trajec- 
tories (and vice versa). In accordance with the charge 
conservation law, at the superconducting region (S) the 
Cooper pairs are formed and the normal current converts 
into the supercurrent. For the perfectly transparent N/S 
interface the charge doubling due to Andreev reflection 
results in enhancement of the subgap conductance by a 
factor of two compared with the corresponding normal 
state conductance. Being the two-particle process An- 
dreev reflection is strongly suppressed due to QP scat- 
tering at the layer of insulator separating the N and S 
regions. Indeed, in case of the small intcrfactial bar- 
rier transparency T C 1, the Andreev reflection proba- 
bility, and therefore the conductance is proportional to 
T 2 which is smaller by the factor T as compared to the 
single- electron case. 

The phenomenon of subgap conductance suppression 
results in the remarkable dependence of transport prop- 
erties of N/S structures on the spatial distribution and 
symmetry of the superconducting gap. For example, re- 
cently the charge current measurements were utilized for 
the direct observation of multi vortex structures in meso- 
scopic superconductors^. Mixed state of mesoscopic su- 
perconductors is formed by a small number of vortices 
and reveals a rich variety of exotic vortex configurations, 
such as vortex molecules and multiquantum giant vor- 
tices, realizing in such samples due to the quantum con- 
finement of the Cooper pairs motion^. In Ref (3 the phase 
transitions between different vortex configurations were 
observed by the multiplc-small-tunncl-junction method, 



in which several small tunnel junctions were attached to 
the mesoscopic superconductor to simultaneously mea- 
sure the charge transport at the different points of the 
sample. The measured transport characteristics were re- 
lated to the local density of states (DOS), depending on 
the local supercurrent density and hence on the configu- 
ration of the vortex system. Generally, when there is a 
uniform supercurrent flowing at the superconductor, the 
excitation spectrum acquires a Doppler shift by the value 
v s pf, where v s is a superfluid velocity, is a Fermi mo- 
mentum. In this case, the minimal excitation energy is 
given by E min = A — v s pf- At large distances from the 
vortex core (r 3> £, where £ is a coherence length) the 
Doppler shift model gives quite a good approximation of 
the quasiparticle (QP) spectrum with v s = H/2mr (sec 
Ref. [J]). As long as the superfluid velocity is small com- 
pared to the critical value Ao/pp, the Doppler shift of 
the gap edge results in the reduction of the height and 
broadening of the superconducting DOS peak^. Close to 
the vortex core (r ~ £) where v s ~ Aq/pfi the gap edge 
E m in goes to zero. However, in this case the Doppler 
shift model does not hold: it completely misses one of 
the remarkable features of the vortex state: the pres- 
ence of low-energy QP states localized within the vortex 
core. These vortex core states were found in the pioneer- 
ing work by Caroli- de Gennes and Matricon (CdGM)£ 
within the more rigorous approach based on the quasi- 
classical solution of Bogoliubov- de Gennes (BdG) equa- 
tions. It was shown that QP states with energy lower 
than the bulk superconducting gap value Ao are local- 
ized within the vortex core and have the discrete spec- 
trum £^ as a function of the quantized (half-integer) an- 
gular momentum fj,. This spectrum of the CdGM states 
varies from Ao to — Ao as (i changes from — oo to +oo, 
crossing zero when [i changes its sign. At small energies 
|e| <C Ao the spectrum is given by e M ps —fxeo, where 
kp = pf/H and eo = Ao/(fei?£). For conventional super- 
conductors, the interlevel spacing e is much less than the 
superconducting gap Ao since (kp£) ^> 1, therefore the 
CdGM spectrum can be considered continuous as a func- 
tion of the impact parameter of quasiclassical trajectory 



2 



b = —[ijkp. The presence of the QP states bounded in 
the vortex core was confirmed in scanning tunnel spec- 
troscopy (STS) experiments by the observation of the 
zero-bias conductance peak at the vortex corei. The 
analogous effect of the zero-bias conductance enhance- 
ment due to the resonant tunneling into the midgap sur- 
face states was thoroughly studied in high temperature 
D-wave superconductors^. 

The goal of our work is to develop a theory to calculate 
the zero-bias conductance of N/S junction in case when 
the external magnetic field generates vortices in the S 
region near the N/S interface. We consider the charge 
transport across the direction of applied magnetic field. 
The N and S regions are assumed to be separated by 
the interfacial barrier, suppressing the Andrcev reflec- 
tion and the electron transport. We predict the strong 
enhancement of the Andreev reflection due to the reso- 
nant tunneling of QP through the barrier into the CdGM 
states localized within the vortex cores. The essential 
physics of this effect is analogous to the one which takes 
place in double-barrier resonant tunneling diode^. The 
resonant tunneling in double-barrier quantum well struc- 
tures occurs if the energy of incident QP wave coincides 
with the resonant energy, then the reflection probability 
is effectively suppressed due to the interference of the QP 
waves within the quantum well. In our case the quantum 
well is represented by the vortex core and the bounded 
low energy QP state consists of coupled electron and hole 
waves of almost the same amplitude. Therefore, if the 
incident electron has resonant energy, then the reflected 
electron wave is suppressed and the hole wave leaking 
from the vortex core contributes to the Andreev reflection 
probability. The important difference between our situa- 
tion and the double- barrier diode case is that the spec- 
trum of bound states is very dense with the characteris- 
tic intcrlcvcl spacing eo ~ Ao/(fcp£) being much smaller 
than the bulk energy gap Ao. At the same time, the 
broadening of this levels due to the finite barrier trans- 
parency can be estimated as SE ~ A T e~ 2a/? , where T 
is the transparency of the interfacial barrier and the fac- 
tor e~ 2a /^ is due to the exponential decay of subgap QP 
at the superconducting slab of the thickness a, which is in 
fact the distance from the vortex to the N/S interface, £ is 
the superconducting coherence length. Hereafter in this 
paper we will neglect the discreteness of the bound energy 
levels assuming that T e~ 2a ^ 3> (kp^)^ 1 . In fact this 
condition is not very restrictive since fcp£ is large in many 
superconducting materials, for example kp£ ~ 3 • 10 2 in 
Nb and ~ 10 4 in Al. Neglecting the discreteness 
of the spectrum of bounded QP states we can use the 
quasiclassical approximation of QP quantum mechanics 
(see e.g. RefQjj). Within such approximation QP move 
along linear trajectories, i.e. the straight lines along the 
direction of QP momentum n = kpfc^ 1 = (cos 9 P , sin 9 p ). 
Note, that for the N/S point contacts of atomic size in 
magnetic field it is necessary to take into account the 
non-quasiclassical divergence of the electron and hole 
trajectories 1 ^. In the present work we consider the trans- 



port properties of wide N/S junction so that its trans- 
verse dimension L y is much larger than the Fermi wave 
length L y ^> Xp = 2n/kp, therefore the trajectory di- 
vergence can be neglected. Also we assume that L y is 
much larger than the distance from the vortices to the 
N/S interface a and the characteristic intervortex dis- 
tance L v . The important point is that the total QP re- 
flection probabilities can be found as a sum of individual 
vortex contributions. Indeed, the intervortex distance is 
much larger than the Fermi wavelength since L v > £ and 
kp£ ^> 1. Therefore, the quasiclassical trajectories (ex- 
cept those with 9 p = ±tt/2) can pass through not more 
than one vortex core. For some directions of QP mo- 
mentum two resonant trajectories (i.e. passing through 
the vortex cores) are coupled by the normal reflection at 
the interfacial barrier. Assuming the specularly reflect- 
ing barrier this coupling occurs for the momentum angles 
at the narrow angle domains near p = arctan(rtL„/2a), 
where n is integer. The width of the resonant angle do- 
mains 89 ~ Te~ 2a /^ <C 7r is determined by the width of 
the resonant vortex core levels. Therefore, the contribu- 
tion of such trajectories to the amplitude of the reflected 
QP waves is negligible. Then, evaluating the amplitude 
of the QP wave we can separate the resonant trajecto- 
ries coupled with bounded QP states localized within the 
different vortex cores. Since different trajectories do not 
interfere with each other, the resulting reflection prob- 
abilities and the conductance can be found as a sum of 
contributions from individual vortices. 

The dimcnsionlcss conductance (further we will mea- 
sure it in terms of the conductance quantum e 2 /irh) in- 
duced by a single vortex at zero temperature can be es- 
timated as follows: G v ~ N r e~ 2a/ ^T. Here N r ~ k F £ 
is the number of transverse modes of the N/S junction 
which effectively interact with the QP states bounded 
within the vortex core (£ is the characteristic transverse 
size of the vortex core). The factor e~ 2a ^T is the one- 
particle tunneling probability through the barrier con- 
sisting of the insulating layer at the N/S interface and 
the superconducting slab of the thickness a. The total 
vortex-induced conductance is the sum of the individual 
vortex contributions G„t = n v G v , where n v = L y /L v is 
the total number of vortices near the N/S interface, L v 
is the intervortex spacing and L y is the transverse size 
of the junction. The resonant mechanism of Andreev re- 
flection exists along with the usual non-resonant scheme 
involving the two-particle tunneling through the interfa- 
cial barrier with the probability T 2 . The corresponding 
zero-bias conductance can be estimated as Go ~ NqT 2 , 
where -/Vo = kpL y /n is the total number of transverse 
modes in N/S junction. Therefore, the total conductance 
of the N/S junction in magnetic field can be evaluated as 
follows: 

G = aN a T 2 + n v f3N r e~ 2a/ tT, (1) 

where the coefficients a, (3 ~ 1. Then, for the strong bar- 
rier T <C 1, the conductance induced by the vortex chain 
with spacing L v becomes dominant at a < a Cl where the 
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threshold distance a c ~ (g/2)\n.(£/TL v ) can be much 
larger than the vortex core size and the coherence length. 

The parameters of the vo rtex lattice can be estimated 
as follows: a, L v ~ y/fo/B, where B is the average mag- 
netic field of the superconducting sample and <f>o = fthc/e 
is the flux quantum. Then, the magnetic filed depen- 
dence of the vortex-induced second term in Eq.([T]) con- 
ductance is given by: 

G vt ~ N 0] ^-e~ 2 ^^T, (2) 

where H C 2 ~ 4>o/£, 2 is the upper critical magnetic field 
of the superconductor. Note, that at zero magnetic field 
B = the conductance is G ~ NqT 2 and at B = H C 2 
the conductance is much larger: G ~ NqT. Therefore, we 
suggest that if the surface barrier is high T ^> 1, it should 
be possible to observe in experiment the field-induced 
increase of the conductance according to the Eq. © . 

In the present work we do not take into account the 
subgap QP states which exists near the surface of su- 
perconductor due to the Mejssner screening of external 
magnetic field 1 ^. Such QP levels lie higher at the energy 
scale than the CdGM states if the density of supercurrent 
near the surface is less then the critical value. Therefore, 
as long as the zero-bias conductance is considered the 
influence of the surface QP states can be neglected. 

The paper is organized as follows. In Sec. [Til a de- 
scription of the model and the basic equations are given. 
In Sec. Mil we solve the scattering problem to find the 
Andreev and normal reflection probabilities. Sec. IIVI is 
devoted to the conductance calculation and Sec|V]to the 
discussion of obtained results. Finally, conclusions are 
given in Sec. IVI1 

II. MODEL AND BASIC EQUATIONS 

Shown on the Fig.([TJ) is the scheme of the N/S junction 
with vortex lines in the S region parallel to the N/S inter- 
face. For the sake of simplicity we assume that there is 
only one quantized QP mode in z direction and take into 
account only the QP motion in xy plane, perpendicular 
to the vortex lines. 

Considering one vortex from the array, the coordinate 
system is chosen so that the z axis coincides with the 
vortex line and the origin at the xy plane coincides with 
vortex phase singularity point. Neglecting the suppres- 
sion of the superconductivity near the N/S interface due 
to the proximity effect we assume that at x > —a (su- 
perconducting region) the order parameter can be taken 
as follows: 

A(r)=A ^(r)e i *«, (3) 

Here Ao is the gap value far from the vortex core, D v (r) 
and $(r) are the dimensionless profile and the phase of 
the order parameter. The particular form of D v (r) is not 
essential for our consideration, therefore it can be chosen 




FIG. 1: Geometry of N/S junction and sketch of QP trajecto- 
ries. The width of the junction is L y . External magnetic field 
directed along z axis introduces vortices in superconductor. 
The distance from the first vortex chain to the N/S interface 
is a, and the intervortex spacing is L v . The electrons are in- 
jected with the incident angle 8o and experience Andreev and 
normal reflection at the N/S interface. 



similar to the model profile of the isolated vortex core 1 ^: 
D v (r) = r/y/r 2 + £ 2 , where £ is the coherence length. 
The phase distribution $(r) consists of a singular part 
$«(r) = arg(r) and a regular part $ r (r), determined 
by the particular mctastablc vortex lattice configuration 
realizing near the boundary. 

In principle, the method developed in the present pa- 
per is applicable to the arbitrary order parameter phase 
distribution corresponding to the metastable vortex con- 
figuration. At first, we solve a generic problem of the 
influence of a single vortex near the N/S surface on the 
zero-bias conductance of the junction. The vortex stabil- 
ity condition given by the London model, requires van- 
ishing regular part of the superfluid velocity at the vortex 
position: (V<f> r — (2-7r/0 o )A)(r = 0) = 0, where A is the 
vector potential. On the microscopic level, this condition 
is necessary for the existence of the CdGM QP states 
forming the vortex phase singularit y 14 ! 15 . The next step 
is the summation of the individual vortex contributions 
to the conductance which add independently. The par- 
ticular vortex configuration near the boundary depends 
on many factors, such as the random pinning potential, 
geometry of the superconducting sample, magnetization 
history, etc. To estimate the dependence of the conduc- 
tance on the magnetic field we consider the model sit- 
uation assuming that the vortices near the boundary of 
superconductor sit periodically on chain with an intervor- 
tex spacing L v at a distance a from the N/S interface. 
We take a and L v as external parameters of the order 
a, L v ~ yj <f>o/B, where B is the average magnetic field of 
the superconducting sample. The influence of the next 
vortex chains on the conductance can be neglected due 
to the rapid decay of the QP tunneling probability with 
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the growing distance from vortices to the N/S interface. 

The expression for the dimensionless zero- bias conduc- 
tance of the N/S junction can be written as follows(see 
RefQ): 



G= 



No r /2 



tt/2 



(1-R n (9o) + Ra(8o))cos6od6 , (4) 



where R n (9o) and R a (6o) are the probabilities of nor- 
mal and Andreev reflection respectively, 0q is the inci- 
dent angle: kp = kp{cos9o, sin#o)- The total number 
of propagating modes N is determined by the channel 
width: N = kpL y /ir. The problem of QP scattering at 
the N/S interface is formulated within the BdG theory. 
The equation for the electron and hole waves coupled by 
the superconducting gap A(r) reads as follows: 




A(r) 
-ffo(r) 



*(r) = £#(r), 



(5) 



Here 



H (r) = — p--A -E F 
2m \ c 



V{x) 



with p = -iW, ^(r) = (u(r),v(r)). The intcrfacial bar- 
rier separating the N and S regions is modeled by the 
repulsive delta potential V(x) — H8{x) 1 parameterized 
by the dimensionless barrier strength Z = H/hVp^. The 
boundary conditions at the N/S interface are: 



*(-a) =0, 



a a *(-o) =(2k F Z)4r(-a), 



(0) 



(7) 



where [f(x)\ = f(x + 0) - f(x - 0). 

To overcome the complexity of the scattering prob- 
lem coming from the broken spatial invariance of the 
superconducting gap, we treat the Eq.© within qua- 
siclassical approximation. Generally, the quasiclassical 
form of the wave function can be constructed as follows: 
#(r) = e ikF - r /(r), where /(r) = (U(r),V(r)) is a slow 
varying envelope function. Then the system ([5]) reduces 
to the system of the first-order quasiclassical equations 
along the linear trajectories defined by the direction of 
the QP momentum n = kpfc^ 1 = (cos 9 P , sin 9 P ). Each 
trajectory is specified by the angle 6 P and the impact 
parameter b = rsin((? — 9 P ), where 9 is the polar angle: 
r = — r(cos 9, sin 9). Introducing the coordinate along tra- 
jectory s = (n-r) = —rcos(8 p — 9) we arrive at the follow- 
ing form of the quasiclassical equation: H f = ef, with 
the hamiltonian: 



where e = E/Aq, £ = fiV^/Ao is the coherence length, 
D v (r) and <E>(r) are the dimensionless magnitude and 
phase of the order parameter and F(r) = (ir£/(j)o)n ■ A. 

Considering the zero-bias problem we will have to ana- 
lyze only the zero-energy excitations with e — 0. Within 
the normal metal region at x < —a, neglecting the in- 
fluence of magnetic field on the QP motion, the quasi- 
classical equation ([5]) becomes trivial: d s f(s,b) = 0. It 
means that the envelope function is constant along the 
trajectories. 

Obviously, this is not the case at the superconducting 
region x > —a, where the electron and hole waves are 
coupled. Note that for the wave functions at the S re- 
gion corresponding to the zero energy the following repre- 
sentation can be used^: / = e c (e^*)/ 2 , e -*('H-<&)/2) ; 
where £ = £(s, b) and r\ = r](s, b) arc real-valued func- 
tions. Then, the quasiclassical equation © can be writ- 
ten as follows: 



£d 8 r] + 2D V cos rj + e d = 0, 



t,d s ( + 2D v sin 77 = 0. 



(9) 



(10) 



where £d(r) = McfVs/Ao is the dimensionless Dopplcr 
shift energy. For the wave functions f± decaying at the 
different ends of the trajectory f±(s = ±oo) = from 
Eci. lfTO"!) we obtain: 



r]±(s = ±oo) = ±7r/2. 



(11) 



As we will see below, the main contribution to the 
enhanced Andreev reflection probability comes from the 
trajectories with small impact parameters \b\ <C £, pass- 
ing through the vortex core. For such trajectories ne- 
glecting the nonsingular part of the superfluid velocity 
near the vortex core the analytical solution of Eq. © can 
be obtained following the results of Reffl6l: 



tan ??± = 0.5 [A^e- 2K ^ ~ A ± e 2K ^) , (12) 
A± = j(b)(sgn(s) =F 1), where 7(6) = — u>b, 



1-1, 



K(s) = jj D v (s')ds' = V(*/0 2 + 



u=- [°° Me-^Wd,. 



III. SCATTERING PROBLEM: NORMAL AND 
ANDREEV REFLECTION PROBABILITIES 



H = —i£tJ z d s + F + D v (a x cos $ — a y sin <f>) . 



The boundary conditions (|6I7[) determine the specu- 
(8) larly reflecting N/S interface, coupling the waves with 
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wave vectors kp = kp (cos 80, sm0o), and kV = 
kp (cos(-7r — o ),sin(7r — 8 )). Therefore if the incident 
electron wave is u, = e lkpr , then reflected electron u r 
and hole v r waves will have the form 



U r e 



V r e 



where and U r (b,s) and V r (s,b) are the envelope func- 
tions. Each point (—a, y) at the N/S interface lies on the 
intersection of two quasiclassical trajectories, character- 
ized by the angles 8 p = 8$ and 8 p = tt — Sq. From the sim- 
ple trigonometry it is easy to see that the impact param- 
eters of these trajectories are 6 + = — asin(#o — 8)/ cos8 
and 6_ = — asin(#o + 8)/ cos 8 correspondingly, where 
= — arctan(y/a) is the polar angle. The coordinate of 
the intersection point is s + = —acos(9 — 80)/ cos 8 and 
s_ = a cos(8 + 8o)/ cos8 for the trajectories characterized 
by the angles 8 P = 80 and 8 P = it — 80 correspondingly. 
Then, the reflection probabilities are given by: 



where x(8, 80) = V- ~ V+ an d Z = Zj cos^ - 

Note, that for the small impact parameters \b±\ <C £ 
the factors e tT,± can be obtained analytically with the 
help of Eq.fTTJ]) as follows: 



, j-j(e-e Q ) 
'j + i(e-8 \ 



.J-i(8 + 8 ) 
J + iid + So)' 



(17) 



(18) 



where J = e - 2K( - a / cos)}S > cos8/(2u>a). For the small an- 
gles \8\, \8o\ <C £/a the equations (jTTJ) and (fT8|) are valid 
simultaneously yielding: 



J 2 + 8 2 - 81 - 2i8 J 
' J2 + 02 _ 2 + 2id J 



(19) 



Rri 



a 



\U r (8,8 )\ 2 (cos 8)~ 2 d8 (13) 



IV. VORTEX-INDUCED ZERO-BIAS 
CONDUCTANCE 



R a (8o) 



L 



\V r (8,8 )\ 2 (cos8)- 2 d8, 



(14) 



y J -a 



where the integration is done over the N/S interface, a = 
arctan(L y /2a), U r (8,8 ) = U r (b-,S-) and V r (8,8 ) = 
V r (b + ,s+). 

Following the usual procedure, to find the reflected 
wave amplitudes U r (8,8o) and V r (8, 9q) one needs to 
match the N and S regions solutions at the N/S interface. 
For the envelope functions the boundary conditions (|6I7|) 
yield: 

1 + U r = e iri +' 2 C + + e ir <-/ 2 C-, 
V r = e- ir >+ /2 C+ + e - lT >- /2 C-, 
(1 - U r ) + 2*Z(1 + U r ) = e" 1+/2 C + - e"'- /2 C, 

V r (l + HZ) = e~"^/ 2 C + - e-' iT '- /2 C~, 

C + ,C _ are arbitrary constants and r]± = rj±(s±,b±), 
where r]±(s, b) are the solutions of Eq.([9]) with boundary 
conditions (fTl"]) along the trajectories with 8 P = 80 and 
Op = 7T— 80 for the upper and lower signs correspondingly. 
Solving this system we obtain: 



U r (8,8 ) 



V r (8,8 ) 



(1 - e^)(Z 2 - iZ) 
l + Z 2 (l-e^) '' 



-«;+ 



l + Z 2 (l-e l x)' 



(15) 



(16) 



Now, using the expressions (fT5l [T6|) for the amplitudes 
of reflected waves and reflection probabilities (fT3l [T4| it 
is possible to find the zero-bias conductance. Introducing 
the function g(8,8 ) = 1 - \U r (8 ,8)\ 2 + \V r (8 , 8)\ 2 the 
expression for the dimensionless conductance ([4]) reads 
as follows: 



h rt f a f IT /2 

G=-£— (cos 8)- 2 d8 g (8, 8 ) cos 8 d8 , (20) 



2 71 



ir/2 



where a = arctan(Ly/2a). It is convenient also to intro- 
duce here a local conductivity, i.e. the conductance per 
unit length of the N/S surface: 

<r(P) = 7T I ' ff(Mo)cos0 o d0 o . 



2- 



tt/2 



Employing Eqs. (|15|16p we obtain: 

g(0,8 ) = 



(Z A + Z 2 )\\ - e^\ 2 + 1 



(21) 



If the applied magnetic field is zero and the super- 
conductor is homogeneous \ — ""j we obtain g(8, 80) = 
go(8 ), where g (8 ) = (l/2)(Z 2 + 1/2)" 2 . Then, the 
vortex-induced part of the conductivity is given by 

kp r /2 

°v(8) = —\ g v (8, 8 ) cos 8 d8 , 



2- 



-tt/2 



where g v = g - g 



9v 



(Z 4 + Z 2 )(4-\l-e^\ 2 ) 
2(Z 2 + 1/2) 2 ((Z 4 + Z 2 )\l - e^p + 1) 
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To start the analysis of Eq. pTj) let's note that for the 
low surface barrier Z — > we get g v (9 7 9 a ) = 0. In this 
case all the incident QP undergo Andreev reflection and 
the zero-bias conductance is the same as in case of ho- 
mogeneous superconductor: G = 2Nq. As the barrier 
becomes higher, the Andreev reflection is suppressed and 
the conductance is reduced. 

The function g v (6,9o) reaches its maximum 

_ o Z^ + Z 2 
9vm ~ (Z 2 + 1/2) 2 

if |1 — e tx \ =0. In fact this condition determines the res- 
onant trajectories, corresponding to the zero-energy vor- 
tex core states modified by the normal reflection from the 
interfacial barrier. The resonant trajectories should pass 
through the vortex core therefore the function g v (8,9 ) 
has a sharp peak at 9o ~ ±#- The width of this peak 
is determined by the barrier strength and the distance 
from the vortex to the surface. For the small angles 
|0o| < C/a with the help of Eq.^ we obtain: 



Within such model the condition of vanishing current 
through the N/S interface (5 X 4> — {2-k / <p )A x ) = is sat- 
isfied automatically, and the vortex stability is achieved 
by setting a = yJfo/B. The numerical plot of the func- 
tion a v (9)/a V Q at different distances a from the vortex 
to the interface is presented on the inset at Fig. ((2]). The 
maximum value a v (9 = 0) with good accuracy coincides 
with the analytical estimation given by (|2"3")1 . The coef- 
ficient in Eq. (|24]) is found to be nearly constant as a 
function of a: it decreases slightly from ~ 0.6 at a = 2£ 
to ~ 0.4 at a = 5£. On the Fig.© for the various mag- 
nitude of barrier strength we plot the ratio a v /ao of the 
average vortex-induced conductivity a v = G v /£ to the 
conductivity of the N/S junction in the absence of vor- 
tices given by: 

k F r /2 

ao = —/ g (9, 9 ) cos 9 d9 Q . 

27T J-n/2 

V. DISCUSSION 



9v(o ,e) 



(Z 2 + l/2)*(0§ -9 2 - j 2 ) 2 + j 2 X 



(22) 



where J = e - 2K ^ / '(2wa) ~ (£/a)e~ 2a /«, which is a 
small parameter since J -C 1 for a > £. The maximum of 
g v (9, 9 ) determined by Eq.® lies at 9 2 = 9 2 + J 2 . 

Employing Eq. (|2"2"|) it is easy to compute the vortex- 
induced part of the conductivity a v (9) at the small angle 
domain \9\ <C £/a. The main contribution to integral 
over 9q comes from the small vicinity of the point 9q = 9. 
Then with good accuracy we obtain: a v = a v o, where 



&v0 — k F Jo 



(Z 2 + 1/2)2 



(23) 



At larger angles the function a v (6) can be evaluated 
only numerically. Numerical calculation described be- 
low shows that u v {9) is maximal at 9 = and steadily 
decreases to zero as \9\ — > tt/2. (see inset on FigEJ . 
Then, the resulting conductance induced by a single vor- 



tex G v = a J" a (cos9) 2 a v (9)d9 is given by: 



G v = 0{k F (,)e~ 2K ^ 



Z 2 



(Z 2 + 1/2) 3 



(24) 



where = (2w£)- x J"Jcos9)" 2 a v {9)/a v0 d9 ~ 1. 

To evaluate the conductance rigorously, we find the 
factor e %x in Eq. (|2"Tj) and then the reflection probabili- 
ties solving numerically Eq. ([9]) with boundary conditions 
(fTTjl . We assume that the regular part of the phase dis- 
tribution is $ r (r) = —arg(r — r av ) corresponding to the 
image vortex situated at the point r av = (—2a, 0,0) be- 
hind the N/S interface. The vector potential is chosen 
as A = B [z x (r — r )] /2, where r = (—a, 0, 0) is the 
point at the boundary between the vortex and the im- 
age vortex, and zq is the unit vector along the z axis. 



Considering the strong barriers Z ^> 1 Eq. (|24j) can be 
written as G v = 0(k F £,) e - 2K< - a ^T 7 where we introduced 
the barrier transparency T = (1 + Z 2 )^ 1 m Z~ 2 . A sim- 
ple understanding of this result can be obtained within 
the framework of the tunneling hamiltonian approach^. 
The conventional expression for the tunneling conductiv- 
ity of the N/S junction at zero temperature reads: 



(25) 



where v is a local superconducting DOS at the Fermi 
level and a n is a normal state tunneling conductivity. 
For the S-wave superconductor the transformation of vor- 
tex core states near the surface can be neglected at first 
approximation^. Then the local DOS near the surface 
at the Fermi level is determined by the density of vortex 
core states v v (r) given b y 19 ' 20 : 



1 

2^ 



\f(r,9 p )\ 2 S(e k F rsin(8 - 9 p ))d9 p , (26) 



where f(r,9 p ) is the envelope of the QP wave function. 
For the CdGM wave functions |/(r, 6> p )| 2 - e~ 2K(r ^k F /C, 
then v v (r) ~ Vo{£, / r)e~ 2K where Vq = m/fi 2 is the 
two-dimensional DOS at the normal metal. Substitut- 
ing v = v v (r = a/ cos 9) and a n ~ Tk F we obtain the 
conductivity: a ~ k F TJ(9), which coincides to the or- 
der of magnitude with expression (|23[) if Z > 1 and 
\9\ <C 1. Integrating a(9) over the N/S interface we ar- 
rive at expression (j2~4")) for the conductance with factor 
given by: 0~ J" a exp(2K (a) — 2K(a/ cos #))(cos 9)~ 1 d9. 
Note, that although yielding the qualitatively right an- 
swer for the vortrex-induced conductance, the tunneling 
hamiltonian approach drops out the contribution from 
the non-resonant Andreev reflection with the probability 
of the order T 2 . As wc will sec below, the contributions 
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of the resonant and non-resonant Andreev reflections can 
be comparable even when the surface barrier is rather 
strong. Certainly, the tunneling hamiltonian approach 
fails to provide the answer if the barrier strength is not 
very high, when the influence of vortices on the conduc- 
tance is reduced and the non-resonant Andreev reflection 
prevails. On the Fig(3^ we show in logarithmic scale the 
vortex induced conductivity function of the bar- 

rier strength Z for the several values of the distance o. 
At small values of Z the function a v (Z) grows a v ~ Z 2 
in accordance with the estimation (f2"5| . At larger values 
of the barrier strength Z ^> 1 the behaviour of a v (Z) 
changes to a v ~ Z~ 2 . But at the same time the conduc- 
tivity without vortices at 2 > 1 behaves as <7q ~ Z~ 4 . 
Therefore, the ratio a v /ao is monotonically growing as 
a function of the barrier strength Z proportional to Z 2 
(see FigHj). 

Let's have a look at the expression for the total con- 
ductance of the N/S junction, which has quite a simple 
form if Z ^> 1. Neglecting edge effects end summing up 
the individual vortex contributions we obtain: 

G = (8/l5)N T 2 + n v (3(k F ^)e- 2K{a) T, (27) 

where n v = L y /L v is the total number of vortices near 
the N/S interface. The obtained expression for the total 
conductance |27|) consists of two terms. The first term 
Go ~ NqT 2 coincides with the conductance of the N/S 
junction at zero magnetic field. The factor T 2 is deter- 
mined by the probability of the sequential tunneling of 
the incident and reflected QP through the high interfa- 
cial barrier. The second term is the total vortex-induced 
conductance G vt = n v G v ~ n v (kp^)e~ 2K ^T; it comes 
from the tunneling of the incident QP into the zero en- 
ergy CdGM states inside the vortex cor e. Indeed it is easy 
to see that e~ 2K{ ^T, where K(a) = ^{a/Cl 2 + 1 - 1, is 
the one- particle tunneling probability through the in- 
tcrfacial barrier and the superconducting layer of the 
thickness a with slightly suppressed gap due to the pres- 
ence of the vortex. The factor kpS, is the number of 
resonant transverse modes for a single vortex. The 
vortex-induced conductance G v t prevails over Go when 
a < a c , where the critical distance a c is determined by 
a c sa (£/2) ln(L v /T£). The parameters of the vortex con- 
figuration, such as the intervortex spacing L v and the 
distance a from the vortex array to the boundary of 
superconductor are determined by the magnetic field, 
therefore the conductance of the N/S junction can be 
controlled by the magnetic field. Using Eq.([2]) we ob- 
tain that the critical magnetic field H c when G v t ~ Go 
is determined by the following transcendental equation 
ln(x/T) = 2/x, where x = yjB/H c2 . Taking for ex- 
ample the barrier strength Z = 5, we obtain that the 
critical field is B c ~ 0.5H C 2, and the critical distance 
a c ~ 1.5£. Therefore, the influence of the resonant vor- 
tex core states on the conductance can become significant 
when the magnetic field is less then the upper critical and 
vortices are quite far from the N/S interface. 

Finally, we should note that in real N/S junctions the 



motion of QP is certainly affected by impurity scattering. 
The influence of impurities can be neglected completely 
assuming that the life time of the vortex core states due 
to the finite barrier transparency is much shorter than 
the relaxation time of the QP momentum: h/SE -C t, or 

l e » (Te- 2 ^)-^, (28) 

where l e = VpT is an elastic mean free path of QP at the 
S region. This condition certainly can bc fulfilled if the 
barrier transparency is not very high, i.e. T ~ 1 and the 
vortex chain is situated not far from the N/S interface, 
so that the factor Te~ 2a ^ is not very small. Otherwise, 
if the condition l]28p is not fulfilled, the impurity scatter- 
ing will modify the conductance. The simplest approach 
to estimate the conductance in this situation is based 
on the tunneling hamiltonian, yielding the Eq. H25|) for 
the conductivity. Due to the impurity scattering the lo- 
cal superconducting DOS differs from that given by the 
Eg. (f2"S|) , which is valid for the clean case l e 3> £■ In par- 
ticular, the sharp peak at r = is smeared, therefore 
at the center of vortex the DOS is smaller as compared 
to the clean case. But at the distances r > £ from the 
vortex core (e.g. at the N/S interface) the DOS is not 
suppressed by the impurities even if l e ~ £. On the 
contrary, it is even larger than in the clean case due to 
the smearing of the DOS peak at the center of vorte:* 2 ^. 
Therefore, in case of the rather high impurity concentra- 
tion on the S side the dependence of the vortex -induced 
conductance on the magnetic field is still described by 
the Eq.([2]). Another important point is the influence of 
impurities on the non-resonant part of the conductance, 
i.e. the first term in Eq.([T]). In particular, the interfer- 
ence of QP waves reflected from the interface barrier and 
impurities on the N side of the N/S junction can also 
result in the low-bias conductance enhancement, known 
as reflectionless tunneling (see Ref. and references 
therein). In experiments where the reflectionless tun- 
neling effect was observed 2 ^ the condition l e > £ was 
fulfilled. In this case the critical value of magnetic field 
suppressing the reflectionless tunneling 2 ^ H c ~ (f) /(12l 2 ) 
was much less than the upper critical field of the super- 
conductor. In our case the applied magnetic field must 
be strong enough to create the dense vortex lattice in 
superconductor: B ~ H c2 » H c . Therefore, under the 
same experimental conditions as in ReflU the reflection- 
less tunneling effect is absent in the range of magnetic 
fields that we arc interested in. 



VI. CONCLUSION 

To summarize, we have investigated the low-energy 
charge transport in N/S junction across the direction of 
applied magnetic field. We have found the strong en- 
hancement of the zero-bias conductance due to the res- 
onant tunneling of the incident QP into the subgap vor- 
tex core states. The effect is most sound if the conven- 
tional channel of Andreev reflection is suppressed by the 
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FIG. 2: Plot of the ratio a v /cro of the average vortex induced 
conductivity to the conductivity in the absence of vortices. 
Curves from top to bottom correspond to Z = 5, 4, 3, 2. Inset: 
function cr v (9)/a v o for a/£ = 2 (open circles), a/£ = 3 (filled 
circles) a/£ = 5 (asteriskes) ; Z = 2. 



high intcrfactial barrier. Note that usually, the vortex 
core states are investigated in STS experiments, where 
the charge transport is measured along the direction of 
magnetic field. For the alternative to the STS methods 
now we can suggest the transport measurements in pla- 
nar structure with wide-area N/S contacts. The vortex 
induced conductance that we have studied depends ex- 
ponentially on the distance from the vortex chain to the 
N/S interface and therefore can be effectively controlled 
by the external magnetic field. Also for the possible ex- 
perimental setup one can consider the mesoscopic super- 
conducting sample with lateral tunneling junctions, such 
as in Ref|2|. Since the vortex -induced conductance is 
proportional to the number of vortices, the conductance 
vs magnetic field dependence should reveal pronounced 
steps marking the switch of the total vorticity of the sam- 
ple. 
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